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We study the intensity distribution function, P(I), for monochromatic waves propagating in 
quasi one-dimensional disordered medium, assuming that a point source and a point detector are 
embedded in the bulk of the medium. We find deviations from the Rayleigh statistics at moderately 
large I and a logarithmically-normal asymptotic behavior of P(I). When the radiation source and 
the detector are located close to the opposite edges of the sample (on a distance much less then the 
sample length), an intermediate regime with a stretched-exponential behavior of P(I) emerges. 
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1. When a wave propagates through a random medium, 
it undergoes multiple scattering from inhomogeneities. 
The scattered intensity pattern (speckle pattern) is 
highly irregular and should be described in statistical 
terms. One of its characteristics is the intensity distribu- 
tion, P(I), at some point r. Almost a century ago Lord 
Rayleigh, using simple statistical arguments, proposed a 
distribution which bears his name: 



P (I) = exp(-I) , 



(1) 



where I is the intensity normalized to its average value, 
I = 1/(1)- The Rayleigh distribution has moments 
(J™) = n\ and it provides, in many cases, a rather accu- 
rate fit to experimental data, as long as / is not too large 
II . For large J, however, the data show large deviations 
from Eq. (H) |l]-|j|. Various extensions of Eq. ([!]) have 
been proposed in the literature. Jakeman and Pusey 
proposed to fit the data with the K-distribution. It con- 
tains a phenomenological parameter rj and its moments 
are given by (I n ) = n\ri~ n T(n + rj)/T(ri). The experi- 
mentally relevant situation corresponds to rj 3> 1. In this 
case moments up to n < rj can be approximated as 



(J n ) ~ nl exp(n 2 /?7) 



(2) 



where only the leading term in the exponent has been 
kept. Thus, only low moments (n -C yjrf) are close to the 
Rayleigh value n\. Some theoretical support to the phe- 
nomenological formula Eq. (|2|) has been given by Dashen 
pj who considered smooth disorder (typical size of inho- 
mogeneities much larger than the wavelength). 

More recently, there was a considerable theoretical ac- 
tivity in studying the statistics of the transmission coef- 
ficients t a b of a one-dimensional sample with short-range 
disorder In this formulation of the problem, a 

source and a detector of the radiation are located out- 
side the sample. The source produces a plane wave in- 
jected in an incoming channel a, and the intensity in an 



outgoing channel b is measured. It was shown in 
that the distribution of the normalized transmission co- 
efficients s a b = t a b/ (t ao ) crosses over from the Rayleigh 
distribution P(s a b) = e _Sab to a stretched-exponential 
one P(sab) ~ e~ 2 ^ gSab , where g is the dimensionless con- 
ductance. 

In this paper we consider a different situation, when 
both the source and the detector of the radiation are em- 
bedded into the bulk of the sample, and calculate the 
intensity distribution P(I) in this case. We prove that, 
for not too large n, the moments can be indeed described 
by Eq.(H), and compute the parameter rj phenomenolog- 
ically introduced in jj]. We compute further the whole 
distribution function P(I) and show that its asymptotic 
behavior at large I is of a logarithmically-normal form, 
in contrast to the stretched-exponential asymptotics of 
P{t ab ) found in Refs. @-||. Finally, we discuss how these 
two different forms of the asymptotic behavior match 
each other and describe physical mechanisms governing 
both of them. 

2. We assume a quasi one-dimensional geometry, i.e. 
consider a tube of transverse dimension W and length 
L » W, filled with scattering medium (Fig. Q). The 
monochromatic source of radiation is placed at point 
r D . The field at some point r is given by the (retarded) 
Green's function Gn(r,r a ) and the radiation intensity is 
defined as I(r, r„) = |G^(r,r )| 2 . The average inten- 
sity (7(r, r )) is represented diagrammatically in Fig |^. 
It consists of a diffusion ladder (a diffuson), T(ri,r2), 
attached to two external vertices. The vertices are short- 
range objects and can be approximated by a (S-function 
times (l/An), so that (J(r,r )) = (£/47r) 2 T(r, r ). For 
the quasi one-dimensional geometry, the expression for 
the diffuson reads 



T(r,r ) = 



3 [ Z< {L- Z> )] 



4tt 



AIL 



(3) 



where I is the elastic mean free path, A is the cross- 



1 



section of the tube, z-axis is directed along the sample, 
z< = min(z, z a ) and z> — max(z,z ). We assume, of 
course, that \z — z | ^> W. 

The intensity distribution P(I), in the diagrammatic 
approach, is obtained by calculating the moments (I n ) 
of the intensity. In the leading approximation [jlO| , one 
should draw n retarded and n advanced Green's func- 
tions and insert ladders between pairs {Gr,Ga} in all 
possible ways. This leads to (I n ) = n\(I) n and, thus, to 
Eq. ©. 

Corrections to the Rayleigh result come from diagrams 
with intersecting ladders, which describe interaction be- 
tween diffusons. The leading correction is due to pairwise 
interactions. The diagram in Fig. || represents a pair of 
"colliding" diffusons. The algebraic expression for this 
diagram is 

C(r,r.) = 2(^)7[n/'.) 

x T(r,n)T(r,r 2 )T(r3,r )T(r 4 ,r ) 

^)| d 3 p[(Vi+V2) . (V3 + V4) 

+ 2 (Vi • V 2 ) + 2(V 3 • V 4 )] f[ S(p -rA , (4) 

»=i J 

where k is the wave number and V 2 acts on r,. The 
factor (^/47r) 4 comes from the 4 external vertices of the 
diagram, the T's represent the two incoming and two 
outgoing diffusons and the expression in the curly brack- 
ets corresponds to the internal (interaction) vertex [ pT| . 
Finally, the factor 2 accounts for the two possibilities of 
inserting a pair of ladders between the outgoing Green's 
functions. Integrating by parts and employing the quasi 
one-dimensional geometry of the problem, we obtain (for 
z Q < z): 



C{z,z )~2(I(z,z )) 2 1 + 



3 n 



(5) 



where (I(z, z )) = z °^ tL ^ is the average intensity, 



7 = 2g 



L 2 (3z + z ) - 2Lz(z + z ) + 2z 2 (z - z ) 



» 1 , 



(6) 



and g = kg£A/3xL 3> 1 is the dimensionless conduc- 
tance of the tube. For simplicity, we will assume that 
the source and the detector are located relatively close 
to each other, so that \z — z a \ <C L, in which case Eq.(|^) 
reduces to 7 = gL 2 /2z(L — z). (All the results are 
found to be qualitatively the same in the generic situ- 
ation zq ^ z — zq ^ L — z ^ L.) 

In order to calculate (I n ) one has to compute a combi- 
natorial factor which counts the number Ni of diagrams 



with i pairs of interacting diffusons. This number is 
Ni = (n!) 2 /[2 2l i!(n - 2i)l] ~ (n!/i!)(n/2) 2i , so that 



(iy 



[ ^ ] 1 /2n 2 



(7) 



Although i cannot exceed n/2, the sum in Eq. (Q) can be 
extended to 00, if the value of n is restricted by the con- 
dition n <C 7. Eq. (Q) represents the leading exponential 
correction to the Rayleigh distribution. Let us discuss 
now the effect of higher order "interactions" of diffusons. 
Diagrams with 3 intersecting diffusons will contribute a 
correction of n 3 /j 2 in the exponent of Eq. (Q), which is 
small compared to the leading correction in the whole re- 
gion 11C7, but becomes larger than unity for n ^ 7 2 / 3 . 
Likewise, diagrams with 4 intersecting diffusons produce 
a n 4 /7 3 correction, etc. Restoring the distribution P(I), 
we find 



P(I)~expf-J+|-7 2 




(8) 



It should be realized that Eq. (g) is applicable only for 
/<7~g and, thus, does not determine the far asymp- 
totics of P(I). The latter is unaccessible by the pertur- 
bative diagram technique and is handled below by the 
supersymmctry method. 

3. In the supersymmetry formalism, averaging over disor- 
der is replaced by functional integration over supermatrix 
fields, Q(r), which satisfy the constraint Q 2 — I p2[ . For 
technical simplicity, we will assume that the time reversal 
symmetry is broken by some magnetooptical effects. The 
integration is done with a weight function exp[— S{Q}], 
where S{Q} is the cr-model action, 



S{Q} 



■nvD 



d 3 rStr(VQ) 5 



(9) 



Str denotes the supertrace, D is the diffusion constant, 
and v is the average density of states. In the consid- 
ered quasi-lD geometry, ttvD = gL/2A, and the field 
Q depends on the z-coordinate only, yielding S{Q} = 
— (gL/8) J dz Str (dQ jdz) 2 . Following the derivation out- 
lined in [pl| Hq| , the moments of the intensity at point r 
due to the source at r Q are given by 



in = 



I6tt 2 



[DQ] (Q b &(z)) n (Q^(z Q )) n e~ s{Q \ 

(10) 



where QiKQai) 1S the retarded-advanced (resp. 
advanced-retarded) matrix element in the boson-boson 
sector of Q. Assuming again that the two points r and 
r Q are sufficiently close to each other, \z — zq\ <C L 
and taking into account slow variation of the Q- 
field along the sample, we can replace the product 



2 



QiaWaaC*)) by Qf 2 {z)Q b ^[{z). We get then the fol- 
lowing result for the distribution of the dimensionless 
intensity y — (167r 2 /k^)!: 



discuss, how these two results match each other. Ana- 



P(y) = J dQ5(y + Q? 2 Ql\)Y(Q) , 



(11) 



where Y(Q) is a function of a single supermatrix Q de- 
fined as follows EMI: 



Y(Q ) = 



Q(r 



[VQ] exp[-S{Q}] 



(12) 



In general, the function Y{Q) depends only on the pa- 
rameters l<Ai<oo,— 1<A 2 <1 entering the stan- 
dard parametrization of the Q-matriccs [pl| . Performing 
the integration over the other degrees of freedom, we find 



d_ 

dy ^ dy 2 
Ai + A 2 



Ai — A2 



d\id\2 
Y(X u \ 2 )6(y + l-Xl) 



(13) 



The evaluation of Y(Q a ) = F(AJ, A 2 ) involves, by its def- 
inition ([l2|), an integration over all supermatrix fields, 
which assume a given value Q a at point z and sat- 
isfy the boundary conditions Q\ z= o t L = -A, where A = 
diag{l, 1, — 1, — 1}. Since g ^> 1, this calculation can 
be done by the saddle point method, as suggested by 
Muzykantskii and Khmelnitskii |T^] , The result is |ll| 



y(A 1 ,A 2 )^exp{-^[^ + 2 2 ]} 



(14) 



where Ai = cosh#i, A2 = cos 82 (0 < 6\ < 00, < O2 < 
7r). In fact, the dependence of Y on 9 2 is not important, 
within the exponential accuracy, because it simply gives 
a prefactor after the integration in Eq. Therefore, 
up to a pre-exponential factor, the distribution function 
P(y) is given by 



P(y)~Y(A 1 = v ^T^,A 2 = l)~exp(- 7 2 /2) , (15) 



where 9\ = ln^l + y + \/y)- Finally, after normalizing 
y to its average value (y) — 2/7, we obtain: 



,. N1 ,<;_2 



In 2 



I + 2I/7+ \/2//7 



(16) 



For J « 7, Eq.(|l6|) reproduces the perturbative expan- 
sion (^|) , while for / S> 7 it implies the log-normal asymp- 
totic behavior of the distribution P(I): 



lnP(/)~-(7/8)ln 2 (8// 7 ) 



(17) 



4. The log-normal "tail" (17) should be contrasted with 
the stretched-exponential asymptotic behavior of the dis- 
tribution of transmission coefficients |t]-|| . Let us briefly 



lyzing the expression for the moments (10), we find that 
when the points z and zq approach the sample edges, 
zq = L — z <C L, an intermediate regime of stretched- 
exponential behavior emerges: 



In P(J) 



-I 
-2\ 



-Lf 2 

3fl 



-l^ln 2 

8z 



16(f) L g 



Kg 

9«I«9^) 2 
(18) 



Thus, when the source and the detector move toward 
the sample edges, the region of validity of the stretched- 
exponential behavior becomes broader, while the log- 
normal "tail" gets pushed further away. In contrast, 
when the source and the detector are located deep in 
the bulk, zq ~ L — z ~ L, the stretched-exponential 
regime disappears, and the Rayleigh distribution crosses 
over directly to the log- normal one at / ~ g. 

Let us now describe the physical mechanisms standing 
behind these different forms of P(I). The Green's func- 
tion G R (r ,r) can be expanded in eigenfunctions of a 
non-Hcrmitean (due to open boundaries) "Hamiltonian" 
as G R (r ,r) = EM T o)M r )(^ ~ E i + hi)' 1 - Since 
the level widths 7, are typically of order of the Thouless 
energy E c ~ D/L 2 , there is typically ~ g levels con- 
tributing appreciably to the sum. In view of the random 
phases of the wave functions, this leads to a Gaussian dis- 
tribution of G fl (ro,r) with zero mean, and thus to the 
Rayleigh distribution of J(r ,r) = \G R (r , r)| 2 , with the 
moments (I n ) = n\. The stretched-exponential behav- 
ior results from the disorder realizations, where one of 
the states ipi has large amplitudes in the both points ro 
and r. Considering both tpi(r ) and tpi(r) as independent 
random variables with Gaussian distribution and taking 
into account that only one (out of g) term contributes 
in this case to the sum for G R , we find (J n ) ~ n\n\/g n , 
corresponding to the above stretched-exponential form of 
P{I)- Finally, the log- normal asymptotic behavior corre- 
sponds to those disorder realizations, where G R is dom- 
inated by an anomalously localized state, which has an 
atypically small width 7; (the same mechanism deter- 
mines the log-normal asymptotics of the distribution of 
local density of states, see Rcfs. |p^| , |l8| ). 
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FIG. 1. Geometry of the problem. Points r = (x ,y ,z ) 
and r = (x,y,z) are the positions of the source and of the 
observation point respectively. 





FIG. 2. Diagram for the average intensity. The diffusion 
ladder is inserted between two solid lines which represent the 
average Green's functions. 
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FIG. 3. Diagram for a pair of interacting diffusons. The 
external vertices contribute the factor (^/47r) 4 . The shaded 
region denotes the internal interaction vertex, see Eq. 
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